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Abstract 

Scattering by (a) a single composite scatterer consisting of a concentric arrangement of an 
outer N-slit rigid cylinder and an inner cylinder which is either rigid or in the form of a thin 
elastic shell and (b) by a finite periodic array of these scatterers in air has been investigated 
analytically and through laboratory experiments. The composite scatterer forms a system of 
coupled resonators and gives rise to multiple low frequency resonances. The corresponding an- 
alytical model employs polar angle dependent boundary conditions on the surface of the N-slit 
cylinder. The solution inside the slits assumes plane waves. It is shown also that in the low- 
frequency range the N-slit rigid cylinder can be replaced by an equivalent fluid layer. Further 
approximations suggest a simple square root dependence of the resonant frequencies on the num- 
ber of slits and this is confirmed by data. The observed resonant phenomena are associated with 
Helmholtz-like behaviour of the resonator for which the radius and width of the openings are 
much smaller than the wavelength. The problem of scattering by a finite periodic array of such 
coupled resonators in air is solved using multiple scattering techniques. The resulting model pre- 
dicts band-gap effects resulting from the resonances of the individual composite scatterers below 
the first Bragg frequency . Predictions and data confirm that use of coupled resonators results in 
substantial insertion loss peaks related to the resonances within the concentric configuration. In 
addition, for both scattering problems experimental data, predictions of the analytical approach 
and predictions of the equivalent fiuid layer approximations are compared in the low- frequency 
interval. 

Keywords Sonic crystals, Scattering, Helmholtz resonance. Elastic shell, Band gaps 

1 Introduction 



Infinite periodic structures do not support wave propagation in certain frequency intervals known 
as stop bands (band gaps) [HE]- For finite arrays of periodically arranged scatterers, referred to 
as Sonic Crystals, the stop bands correspond to the frequency intervals of very low transmission 
[H m IS]. The position of the stop bands can be tuned to selected frequency intervals by changing 
the spacing between the scatterers. This makes Sonic Crystals attractive for applications as noise 
barriers particularly for narrow band sources. 

The performance of the Sonic Crystals can be improved by increasing the filling fraction [21 [6], 
varying the arrangement of the scatterers [7] and also by replacing the scatterers with the resonant 
elements [U El [lOl HH ESI [l3l El [El [E] . The latter allows increasing insertion loss in the low- 
frequency range if the resonant frequencies of the scatterers lie below the first stop band associated 
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with the array periodicity. The most obvious choice for such resonant elements are spht ring 
resonators [U [131 [15] which essentially are 2D Helmholtz resonators. Alternatively an array of thin 
elastic shells in air has been shown to possess strong multiple resonances in the low-frequency range 

mm- 

Arrays of composite scatterers such as concentric split ring resonators [H |T7] are shown to 
support band gaps generated by the multiple resonances of each scatterer in the array. In this 
paper we investigate the concentric arrangement of two types of the resonant elements i.e. thin 
elastic shells and split ring resonators with multiple slits, referred to as N-slit cylinders. These 
composite scatterers, referred to as composites, serve a double purpose. First, they support multiple 
resonances including Helmholtz and annular cavity resonances and also axisymmetric resonance of 
the elastic shells. The latter is defined by the geometry of the shell and not by the area of the cavity 
that therefore gives advantage in tuning of the composite multiple resonant behaviour compared 
to that for the concentric split ring resonators. It is also shown that the coupling between two 
resonators leads to the shift of the axisymmetric shell resonance to the lower frequency range. 
The analogous effect for coupled mechanical oscillator (multiple degree of freedom system) is well 
known [18]. Second, for the purpose of creating a noise barrier an array of elastic shells exposed to 
the environment would not be practical. A simple solution is to protect the shells by using outer 
cylinders containing slits surrounding inner elastic shells arranged concentrically. 

The scattering problem for the hollow cylinder with non-uniform boundaries has been solved 
semi-analytically in |19[I20] and references therein. This involves solving the linear algebraic system 
of equations derived from the analytical evaluation of the integral equations. An alternative ap- 
proach has been developed in grating theory [21[I22] where solution of the studied partial differential 
equations is subject to the boundary conditions dependent on the coordinates. 

Three configurations involving the N-slit cylinder are studied in this paper (a) alone, (b) con- 
centric with an inner rigid cylinder and (c) concentric with an inner elastic shell. To solve the 
scattering problem for N-slit cylinders we use boundary conditions dependent on polar angle |21] . 
In section [2] this approach is employed to solve a single composite scatterer problem that is similar 
to that solved for an electromagnetic case |22j . However in the present model the finite thickness 
of the wall of the N-slit cylinder is taken into account. This avoids the use of adjustable param- 
eters introduced for numerical stability |21[ [22] and provides a more accurate description of the 
real structure used in the experiments reported later. The solution inside slits is replaced by jump 
conditions that describe the slit interface as a moving piston. It should be noted that use of 
jump conditions to replace the slits makes it possible to accurately predict the total wave field in 
the low- and mid-frequency range. 

It is also shown that in the low-frequency range the slit cylinder can be replaced by an equivalent 
fiuid layer. The approximations are based on the results derived for the perforated plate [M l[25ll26] . 
The developed approximations are of particular interest for studying the low-frequency resonant 
behaviour of the single composite arrangement and a finite periodic array of them. This behaviour 
is attributed to the Helmholtz-like resonances with wavelength much bigger than the characteristic 
sizes of the cylinder (i.e. radius and width of the slits). 

Experimental validation of the analytical and numerical predictions has been carried out. In 
particular the results obtained with the grating theory are compared with the low-frequency ap- 
proximations for the cylinders with different number of slits and various inner structure (N-slit 
cylinder alone, rigid core and elastic shell core). However the analytical models could be used also 
to study the acoustical properties of arrays incorporating resonators of various thickness and with 
different slit width. 

In section [3] the model is generalised for finite arrays of N-slit cylinders and composite elements. 
In section H] the experimental setup is described and the model predictions are compared with the 
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(a) 



(b) 



Figure 1: (a) Cross-section of composite element consisting of a concentric arrangement of an outer 
4-slit rigid cylinder and an inner elastic cylindrical shell, (b) Geometry of the slit. 

data. 

2 Single scatterers 
2.1 Formulation 

Consider the two-dimensional problem of acoustic wave scattering by a single N-slit rigid cylindrical 
shell of thickness h and external radius Tq. The sound is generated by the cylindrical point source 
placed at the origin of the coordinate system which is defined by either the Cartesian (x, y) or polar 
(r, 9) coordinates. Throughout the paper the time-harmonic dependence is taken as exp{—iu}t). The 
widths of consecutive slits in Oxy plane are denoted by d„, re = l..iV and they are infinitely long in 
the direction of the cylinder main axis Oz, see Figure [TJl^a) . The external length of the rigid arc is 
given by D„. It is also assumed that the thickness of the rigid shell is much smaller than its radius 
so that the following geometric simplifications can be applied: 

(a) Internal and external arc lengths of the rigid strip have the same length; 

(b) Internal and external arcs subtend the same angle. 

(c) The angle 0„ subtended by the arc of n-th slit and its width dn are related by 

dn = 4>nro, n = 

(d) The angle ■i/'m subtended by the arc between slit (n-1) and slit n and its width Dm 
are related by Dm = i^mTo-, 'm = 1..N 

Without loss of generality it is assumed that radius of the vector passing through the middle of the 
first slit makes either zero or it angle with the Ox axis. 

For simplicity, the acoustic environments outside, inside of the N-slit cylinder and in the slits 
are assumed to be identical and are described by density p and sound speed c. We consider three 
configurations: (a) an N-slit cylinder alone, (b) with a rigid cylindrical core inside and (c) with an 
elastic shell inside. 

The displacement potential p{r) in the acoustic medium satisfies the Helmholtz equation 



Apa + k^Pa = 0, 



(1) 



where k = u/c is the wavenumber, index a relates solution p to one of the regions (i.e. 'o' 
outside the configuration, 'i' inside it and 's' within the slits) and Laplacian A is given by either 
1 d f d\ 1 92 92 92 
r dr \ dr J r"^ 89'^ dx"^ dip' 

The outer solution "Poir^ ^-^so satisfies Somerfeld's radiation conditions 

-J^ — ik-po = o , as r — )• 00, (2) 

where r = sj ^ ip . 

To proceed with the boundary conditions imposed on the surface of the N-slit cylinder we first 
solve an auxiliary problem for the slits. The geometry of this problem is illustrated in Figure [T]^b) . 
Waves propagating inside the slit are described by plane wave solution as 

= Cxf^^y + C2e-'^\ (3) 

where C\ and C2 are unknown coefficients. It is assumed that contribution of the x-dependent 
components to ps is negligible when slits are small compared to the wavelength (that is kdn *C 1 
and kh<^l). 

Solution ([3]) is subject to the continuity conditions imposed on the slit faces y = ±h/2 as 



Pa =P. 

where qa = dpa/dn is the normal derivative on the slit faces and index a = o^i. Using solution ([3]) 
and its unknown coefficients Ci and C2, equations @ can be transformed to the jump conditions 
[23] 

Pi=Po- hqo, 

Qi = k^hpo + go, (5) 

that relate the wave field inside the N-slit cylinder to that of the outer region at the slit interface. 
Note that due to assumption kh <^ 1 the trigonometric functions in ([5]) are replaced by their leading 
orders. 

The solution of equation ([1]) is subject to the jump conditions ([5]) as well as to the Neumann 
conditions imposed on the rigid surface of the N-slit cylinder. The former and the latter can be 
combined into the set of two boundary conditions |21] that are 



or or 

-T. T— iPo - Pi) = 0, (6) 

Or n 

where stepwise function f{6) of the polar angle 9 introduces the distribution of slits along the N-slit 
cylinder surface over the interval 9 G [0, 2tt] as 

f(9) = S^ if [-0i/2,,^i/2]U..U[2^-0i/2-,^jv-V'iV,2^-0i/2-Viv], 

^ ' \ otherwise. ^ ' 

For the non-zero value of this function equations ([6]) are transformed to the jump conditions ([5]) 
whereas the Neumann conditions can be obtained from ([6]) by setting f{9) to zero. 
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The periodic angle distribution of N slits can be introduced through the alternative form of 
f{e) that is 



1 

+ ^ (0 - 27rj + 4>i/2) -H{e- iTTj - 0i/2)} , (8) 

j=0 

where H{6) is the Heaviside function. 

The total wave field outside the slit cylinder is represented by Po which can be found as 

+00 

Po = <)(A:r)+ ^ ^„i7W(fef)e-^ (9) 

n=— 00 

where r(r, 6) and 9{r, 9) are the polar coordinates of receiver with origin placed at the centre of the 
scatterer and An are unknown coefficients. The solution inside the slit cylinder is given by 

+00 

P^= [BnJn{kr)+CnYn{kf)]e''^\ (10) 

n=—oo 

where coefficients Bn and C„ have to be found. The coefficient Cn can be derived in terms of Bn 
or set to zero according to the type of the core that results in the following 

Cn = BnCn, n G Z, (11) 

where 

Cn = 0, N-slit cylinder alone, (12) 

C„ = -|||^l}, rigid c„,.e. (13) 

" " ~ J^{kR)Y^,{kR)[l - {ksRf + n^] + K - {ksRnp{p,,nRh)-^' elastic shell core, (14) 

where primes (') denote the derivative with respect to polar coordinate r. In equation ()14p the 
elastic shell is described by frequency parameter kg = lu/cs, density ps, dilatational wave speed Cs, 
radius of the elastic shell mid-surface R = (ai + a2)/2 and its half thickness hg = (oi — 02) /2. Here 
the dilatational wave speed Cg for a thin elastic plate is given by 



with Young's modulus E and the Poison's ratio 1^. In this paper we consider thin viscoelastic shells 
made of latex [T5] . For the numerical predictions the size of each elastic shell is always defined by 
the outer radius ai = 0.02 m and thickness 2hs = 0.00025 m. 

In order to find the unknown coefficients An we apply Graf's addition theorem [28] to the outer 
solution Q. This enables us to express (r, 6) in terms of (f, 9) [1Q\. Then, substituting ([9]) and ([10]) 
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into the boundary conditions (0) and taking the inner product 



27r 



< • > exp(— 



we arrive 



at infinite algebraic system of equations in An, n ^ 1,, variables as follows 

oo ( 

AJ6„,^n27rhHW\kro) 

n=—oo I, 



(16) 



j=-oo 



i=—oo (_ 

eh 



2tt 



Jn{kro) Y Fj^rnFn-jIj>, m e 



j=-co 



where 5m,n is Kronecker delta, vector Q = Q{cos a,sma) is the radius vector to the centre of the 
slit cylinder, 

Jn{kri) -\- CnXnikri) 



" J'n{kri)+CnY:,{kny 
and Fn is the Fourier transform of function f{9) given by 

f N 



(17) 



for n = 0, 



1=1 

2 sr^ . f ncpi 
— > sm 

1=0 ^ 



(18) 



Note that factor C„ and the geometrical parameters of the slits only appear in /„ and Fn, respec- 
tively. This makes form of the system (|16p invariant with respect to scatterer core type and the 
arrangement and size of the slits. 

Taking (/)/ = 0, / G Z, the solutions of system (|16p is reduced to the case of scattering by rigid 
cylinder that is 



Ar. 



(19) 



To find the numerical solution the infinite system ()16p has to be truncated at some number 
m,n,j = —M..M. In general the convergence of the numerical solution is dependent on the radius 
To of the N-slit cylinder, number of slits and their angles (pi (i.e. their length) as well as on 
the frequency. Reducing the angle (pn of the slit results in faster convergence so that numerical 
solution approaches the value defined in (fT9]) . It is found that for kvo < 5 and (pi < 0.1, I £ TL, 
truncation number M between 30 and 40 gives accurate results to less than three significant figures. 
Throughout this paper we use M = 35. 

2.2 Results 

Knowing the coefficients An, n £7L, enables us to compute scatterer insertion loss given by 



IL = 20 log 



10 



Po 



(20) 
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Figure 2: Semi-analytical predictions ( ) compared with data ( ) for single scatterers. 

Insertion losses are computed at the point (1.5775, 0) m. The rigid shell has outer radius = 0.0275 
m, thickness h = 0.002 m and N identical slits of width dn ~ 0.004 m, n = 1..N,. (a) Single slit rigid 
cylinder (IS); (b) N-slit rigid cylinder with N=4 uniformly distributed slits (4S); (c) Concentric 
arrangement with outer 4S and inner rigid cylinders with ai = 0.011m; (d) Concentric arrangement 
with outer 4S rigid cylinder and inner latex elastic shell. 



with acoustic potential po found by formula Q. 

The insertion loss computed for a single N-slit rigid scatterer is compared with data in Figured 
(the experimental setup is described in section d]). The results illustrated in Figures [2)^a),(b) and (c) 
are obtained for N-slit cylinders alone and with a rigid core. For the given scatterer the prediction 
exhibits Helmholtz resonance in the low-frequency range (i.e. kro < 1). In both semi-analytical 
results and data, this resonant behavior is observed around 600 Hz for a single slit rigid cylinder, 
see Figure [2)^a) , and around 1500 Hz for a 4-slit rigid cylinder, see Figure [2)^b). It is observed that 
increasing the number of slits results in significant shift of this Helmholtz-type resonance toward 
higher frequencies. The resonance can also be shifted by reducing the cavity area. This effect is 
seen in Figure EJc) where the presence of a concentric internal rigid cylinder results in a shift of 
the resonance to a higher frequency / ~ 1700 Hz. 

The second type of the resonant behaviour observed in Figure [2]is related to the cut off frequency 
of the acoustic mode in the circular /annular cavity. For the single slit rigid cylinder this resonance 
is found around 4200 Hz (Figure [2]^a)). The 4-slit rigid cylinder (Figure [2]Jb)) exhibits this type 
of resonance at a slightly higher frequency / ~ 4400 Hz. For the annular cavity formed by the 
concentric 4-slit rigid cylinder and rigid internal cylinder one of the cut off frequencies can be found 
around 3500 Hz (Figure [2jc)). 

Figure [2]^d) illustrates predicted and measured insertion losses for a 4-slit rigid cylinder with 
an inner elastic shell. The presence of the inner elastic shell with radius ao = 0.02 m causes 
resonant behavior to appear around 1000 Hz, 2700 Hz and 3400 Hz. The first two resonances can 
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Figure 3: Solution ([9|)- (|16p ( ) for single N-slit scatterer compared with its approximations 

( ) defined by ([2^ . ([50]) and ([55]) . Insertion loss is computed at the point (1.5775, 0)m. (a) 

IS; (b) 4S; (c) Concentric arrangement with outer 4S and inner rigid cylinders with ai = 0.011m; 
(d) Concentric arrangement with outer 4S rigid cylinder and inner latex elastic shell. 



be related to the interaction of the axisymmetric resonance of the elastic shell and the Helmholtz 
resonance. The frequency of the axisymmetric resonance of the elastic shell, which is originally at 
about 1200 Hz, is decreased to about 1000 Hz by the concentric arrangement. The second resonance 
at / ~ 2700 Hz is related to the Helmholtz resonance. This frequency is higher than that for the 
empty slit cylinder due to the effective reduction in the size of the cavity. 



2.3 Low- frequency approximations 

In the low-frequency range a simpler approximation can be used to describe scattering by N- 
slit cylinders. The approximation is based on the results derived by Horace Lamb |24j for the 
electromagnetic waves transmitted through a metallic grating. Here we only consider a case when 
N identical slits of width dn = d, n = 1..N, are distributed periodically along the scatterer surface. 
In this case N-slit cylinder can be replaced by a homogeneous fluid shell. The latter is characterised 
by two effective parameters that are density pi and speed of sound q . 

Similar approximations have been used for a periodic array of circular scatterers [27j. The 
fluid layer is defined through the thickness hi = hJ- and density pi = p/J-'^ with filling fraction 
J- = Nd/ (27rro) satisfying the following conditions 

and A = 4 > (21) 
4 / 

The parameters hi and pi were obtained by matching transmission coefficients derived for the 
infinite homogeneous and perforated plates. Having thickness as an effective parameter for the 
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infinite plate gives quite accurate results in the low-frequency regime. However the fact that hi is 
not equal to the actual thickness of the element leads to unnecessary complications in the model 
for the cylindrical shell. 

In contrast here the thickness of the equivalent fluid layer is made equal to the actual thickness 
of the N-slit cylinder and, therefore, the inner and outer radii of the corresponding fluid shell 
are uniquely defined through those of the N-slit cylinder. The obtained approximation gives first 
resonance position within 5% of that for the N-slit cylinder provided that length of d satisfies 
conditions (I2ip . 

To find the unknown density pi and speed of sound q we first write down the plane wave 
transmission coefficient of a fluid layer of thickness h as 

(22) 



(1 + Z-^^e-''"'' - (1 - Z-^^e^^i'^ ' pc 

where ki is the wave number in the fluid layer and Z is its relative impedance. Note that coefficient 
T is derived for the plane waves propagating perpendicular to the surface of the layer. In the low 
frequency regime when kh <^ 1 and kih <^ 1 the transmission coefficient can be approximated by 

^ ^ Z-^ - 0.5ikih{l + Z-^)' ^^^^ 

For a perforated plate such that kd <^ 1 the transmission coefficient T is described by oscillating 
piston-type waves propagating in slits. Using existing analytical results [26] we arrive at 

(24) 



(1 + Z'^^e-''^^ - (1 - Z-^ fe'^'' ' 
with 

Z = :F-^^±'^^^ = :F-^kA, (25) 

n=l 

where A is referred to as end correction. The lower limit of the end correction can be estimated 
from [Ml [25] 

A = -log--i— -. (26) 
vr sm(7ry-/2j 

By assuming low frequency and small filling fraction (i.e. ^ 1) the coefficient T is approximated 
as 

T 

T ~ (27) 

T -0.5ikh{l + T^ + 2A/h)' ^ ' 

In order to be able to replace the perforated plate by an effective fluid layer we equate T and 
T. Comparison of these two coefficients gives relative impedance Z and the relationship between 
wave numbers ki and k in the following forms 
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It is then possible to derive the density and sound speed of the effective fluid layer as 

h , , 

Q = 1 tC (29) 

h + 2A 

P, = —p. (30) 

To derive the solution of the appropriate scattering problem for the layered cylinder continuity 
boundary conditions have to be imposed at the faces r = ro,rj of the fluid layer in the following 
form 

P 

Pl = Pa 

Pi 

— = — (31) 
dr dr ' 



where solutions pa are given by (|9]) and pOj) whereas solution inside the fluid layer is defined 
through the new set of constants Dn and En, n G Z, as 

pi= Yl [DnJn{kir)+EMkif)]e'^'. (32) 

n=— oo 

Note that conditions (j3T]) and solution (f32|) replace the boundary conditions dH) imposed at the 
faces of the N-slit cylinder. 

Using boundary conditions (j3T]) the unknown coefficients An, Dn,En and Cn can be found. The 
coefficient An takes the following form 

An = -Z„i/(i)(fc„Q)e— (33) 

where 

^ _ piWiW3J:,{kro) - pJnjkro) [WiJ^^jkirp) + W2K{kiro)] 

within which 

-Hniro) = Hi,^\kro), 

Wi = [-pYn{kiri)Yn{kri) + piYnikiri)Yn{kri)] Cn + piJ'n{kri)Yn{kiri) - pJn{kri)Yn{kiri), 

W2 = [pJnikiri)Yn{kri) - piJn{kiri)Y^{kri)] C„ - piJ'^{kri)Jn{kiri) + pJn{kri)J'^{kiri), 

WiJn{kiro) + W2Yn{kiro) 
^3 , (35) 



with Cn defined in (fT2]) . 

To validate the low-frequency approximation its predictions are compared with those of the 
model presented in section 12.21 In Figure [3] the results are shown for the same types of scat- 
terers. The accuracy of the approximation deteriorates with increase in frequency. However in 
the frequency range that is inside of the limit (|21|) the approximation reproduces the shift of the 
resonances observed for the concentric arrangement in Figures [3)^c) and (d). 

The low- frequency resonance displayed in Figure [3] below 2000 Hz can be found by constructing 
approximation of the appropriate eigenvalue solution. In the eigenvalue problem the outer acoustic 
potential ([9]) is defined through Bessel functions of second kind only. The solution of this problem 
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corresponds to zeros of the determinant of coefficient An in p3]) with T-ln{ro) = Yn{kro). To obtain 
the approximations of the resonant frequencies we need to consider two different cases defined by 
the core type factor Cn in ()12p - (jl4p . For N-sht cyhnders that are either empty or concentric with 
an inner rigid cyhnder the leading order approximations can be derived by assuming that 



kro 
P_ 
Pi 



(36) 



where /3 <C 1 is a small parameter. These four relations define the frequency regime and contrast 
between fluid layer and outer/inner environment. Note that according to equation ()30p and the 
assumption that <C 1, density pi is always much greater than p. 

Expanding Bessel functions and collecting similar orders of smallness the first resonant frequency 
i.e. n = can be approximated by 



/ 




2p 



af) log (ro/rj) 



(37) 



The resulting values of resonant frequency can be further approximated by using the fact that 
ratios h/ri and h/ro are small. This follows from the assumption made in the beginning of section 
12.11 The expansion log(ro/rj) = h/ri + 0{h?' /rf) and substitution of the equation (f30]l give 




Nd 



7r{rf-ai){h + 2A) 



(38) 



For a single slit = 1 equation (|38p reduces to the well-known value of Helmholtz resonance. From 
equation psp it is also seen that increasing the number of slits leads to the shift of the Helmholtz 
resonance toward higher frequencies by approximately \/iV. This is in accordance with the results 
observed in Figures [3] (a) and (b) which show a shift of Helmholtz resonance frequency by a factor 
of 2 as N increases from 1 to 4. 

Next is considered the case of the concentric N-slit cylinder and elastic shell which is described 
by factor C„ in equation (I14p . The interaction of the axisymmetric and Helmholtz resonances leads 
to the considerable shift of the low- frequency resonance as well as the high-frequency resonances. 
The new position of the low-frequency resonance can be approximated by introducing the following 
order of smallness of the physical and geometrical parameters 



kvo 


= 0{/3), 




= 0{fi) 


P_ 

Ps 


= 0{f3') 


Pi_ 

Ps 


= 0{(3) 


Cs 

c 


= om, 


Cs _ 
ci 


--0{(3), 


hs 
R 









(39) 



For n = the leading order approximation of the resonant frequency becomes 



2Tr \ 2hsR{rl - R^) \og{ro/r,) 



P 2 , Cg hg 2 

Ps R 



R' 



log - + 2^hsR 
U pi 



± 



-rf + i^irf - R')) log ^ - 2UsR 

Ps R J n pi 



+ 8^AhsRHog^ 

Pl Ps Ti 



1/2 



(40) 
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In this equation the presence of two sohitions explains the interaction of the first two low-frequency 
resonances. Two limiting cases can be constructed which show the degeneration of the two res- 
onances into single one. It can be seen that if radii and ro tend to infinity the equation (j40p 
transforms to the resonance of the elastic shell under membrane compression loading [16j i.e. 

27r y hs Ps 

Comparison of the resonant frequency predicted by (j4ip for an elastic shell made of latex with the 
smallest value in (j40p obtained for the concentric 4S rigid cylinder and latex elastic shell enables 
us to estimate the shift of the resonance observed in the composite. The numerical value of the 
axisymmetric resonance (j4ip is / ~ 1270 Hz whereas the lowest axisymmetric resonance of the 
composite is / ~ 1060 Hz. Thus a reduction of the frequency of the axisymmetric elastic shell 
resonance by about 200 Hz is predicted. 

Equation (f40]l can also be reduced to the form of the resonant frequency (|37l) by assuming 
higher contrast in ([39|) between the fluid and elastic media. 



3 Arrays of scatterers 
3.1 Formulation 

The method developed in the previous section can now be adapted to finite arrays of scatterers. We 
first formulate the method of multiple scattering that is described by the superposition of scattered 
wave fields of each element of the array [5] and boundary conditions ^ imposed on the surface of 
all scatterers in the array. As a result the solution po of the Helmholtz equation ([T]) in the acoustic 
environment outside scatterers takes the following form 

Po(,r,9) = Hl,'\kr) + Yl E expM„), (42) 

m=l n=—oo 

where M is the number of scatterers in the array, variables fra{r,9), 6mir,0) are the polar coordi- 
nates with origin in the centre of scatterer of index m and are unknown coefficients. The form 
of the solution inside of the m-th slit cylinder is the same as in the case of single scatterer problem 
that is equation pO]) . 

Using Graf's addition theorem we can rewrite solution (j42l) in terms of polar coordinates 
(^m,^m) of m-th scatterer that gives 

Poirm,0m)= Yl <^(^^m)^i'^(fcQm)e-*"("+"-) + A;^F«(fcf^) 

ra=— oo I. 

Af +00 

+ E ^W^^rn)Hi^'hkQrnp)e'^'-''^^''"^'^^ (43) 

where vector Qm = Qm{cos am, sin am) is the radius vector to the centre of m-th scatterer, vector 
Qmp = Qmp (cos amp 1 sin amp) defines the position of p-th scatterer with respect to m-th scatterer 
and Qm, Qmp > that is the requirement of addition theorem. Expansion (I43p along with the 
inner solution (jlOp are subject to the boundary conditions ([6]). This results in the infinite algebraic 
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Figure 4: Semi-analytical predictions (----) compared with data ( ) for 7x3 arrays of 

scatterers. Distances from array to the source and receiver are 1.5 m and 0.05 m, respectively. The 
scatterers are arranged in square lattice with lattice constant L = 0.086 m. (a) Empty 4-slit rigid 
cylinders (4S); (b) Concentric outer 4S cylinder and inner latex elastic shell. 

system of equations in j4™, m = 1..M, n £ 7j, variables given by 

n=—oo I. 



^ n—m 



j=-oo 



oo oo 



s=l,s^pn=—cov=—oo \. 



(44) 



j=-oo 



^ HW(kQp)e-'^^^+''^U -6^,n27ThJ',{kro,p) 

n=—oo I. 

k'^h °° 1 

+ Fn-^ [Jn{kro,p) - Jn{kro,p)I?] - —Jn{kro,p) Yl Pf-mK-j^! ymeZ,p= 1..M, 

j=-co ) 

where factors Im and F^ are given by equations (fT7|) and respectively. The superscript p 

relates these factors to p-th scatterer. To find unknown coefficients we again use truncated system 
of J\f{2M + 1) algebraic equations. Note that for the slit angles (jyf = 0, I £ Z, p = l..J\f, equations 
(j44p are transformed to the well-known algebraic system describing multiple scattering problem for 
an array of rigid cylinders [5]. 

3.2 Results 

Solving the algebraic system of equations (01]) enables us to compute the total wave field (02]) and 
insertion loss ()20p . The latter is analysed in this section. 

Figure Hl^a) shows measured and predicted insertion loss spectra for a 7 x 3 array of 4-slit 
cylinders. The band gap effect related to the cavity resonance (see Figure [2]) and associated with 
a high insertion loss peak is predicted and observed in the vicinity of / = 1500 Hz. This effect is 
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Figure 5: Semi-analytical predictions for 7x3 array of composites ( ) compared with their 

low-frequency approximations ( ). Distances from scatterer to the source and receiver are 1.5 

m and 0.05 m, respectively, (a) 4S; (b) Concentric outer 4S rigid cylinder and inner latex elastic 
shell. 




Figure 6: Band diagram for the infinite periodic array of composites (left inset) compared with 
the low- frequency approximations of their finite 7x3 arrangement (right inset). All geometrical 
parameters are identical to those used in Figure [5j The enclosed shells in the infinite periodic 
structure are taken without viscous effects and their Young's modulus is identical to that used in 
[l6], Figure 6. (a) 4S; (b) Concentric outer 4S rigid cylinder and inner latex elastic shell. 



followed by the insertion loss peak related to the first Bragg band gap which is observed around 
/ PS 2000 Hz. 

In Figure Sl^b) the insertion loss is computed for the array of composite scatterers consisting 
of concentric 4S rigid cylinders and latex elastic shells. The coupled resonances of this composite 
described in section 12.21 generate various band gaps that result in high insertion loss peaks. One 
of these peaks is observed around 1000 Hz and is related to the band gap due to the shifted 
axisymmetric resonance of latex elastic shell. The insertion loss peak associated with the shifted 
Helmholtz resonance is observed around 2600 Hz and it follows the first Bragg band gap peak 
observed around / ~ 2000 Hz. 

3.3 Low-frequency approximations 

Here the model of equivalent fluid layer is applied to the problem of multiple scattering. Using the 
low-frequency approximations (j29l) and (I30p together with the boundary conditions (I3ip one can 
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derive the analogue of the algebraic system of equations (HH) that is 

N oo 
p=l,p^m q=—oo 

where Zq is given by (j34p for the scatterer with index p = 1..N [5]. The truncated linear system of 
N{2M + 1) equations is solved to find the unknown coefficients and, as a result, the insertion 
loss ()20p . The truncation number M is taken between 5 and 7. This results in substantial reduction 
in the computational time compared to that for the full problem (I44p . 

In Figure [5] insertion loss is computed for the same arrays as considered in section I3.2i It is 
observed that replacement of the N-slit rigid cylinder by the equivalent fluid layer gives results 
accurate to within 5% in the frequency interval / G (0, 2000) Hz. Therefore, for the given geometry 
of the array the proposed approximations accurately predict the presence of the band gaps below 
and including the first Bragg band gap. 

The band diagrams plotted in Figure El^a) and (b) are aligned with the insertion loss spectra 
obtained previously in Figure [5l It is demonstrating that the band gaps coincide with the maxima 
of insertion loss. The diagrams also show that the first band gaps (at 1500 Hz in Figure [6^a) and 
at 1000 Hz in Figure E^b)) are complete. This results in angular independent insertion loss peaks 
bounded in the frequency interval of the complete band gaps. 

4 Laboratory measurements 
4.1 Cylinder constructions 




(a) (b) 

Figure 7: (a) Components of the resonating scatterer: (on the left) photograph of PVC pipe showing 
the discontinuous form of a (single) slit and (on the right) photograph of latex elastic shell, (b) 
Cross section of the concentric arrangement of inner latex cylinder and outer PVC cylinder with 4 
symmetrically-placed slits 

2 m long 0.25 mm thick Latex sheets have been formed into cylinders with outer diameter of 
43 mm. This has been achieved by overlapping edges by a few mm and gluing them together. Slits 
have been gouged from the walls of 55 mm outer diameter PVC pipes. To ensure the structural 
integrity of the pipes, rather than making continuous slits along the complete length of the pipes, 
the slits were made in sections of approximately 20 cm length separated by about 2 cm (see 
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Figure [7Ja)). Concentric arrangements of pairs of latex and 4-slit PVC cylinders were formed as 
shown in Figure [7)^b). To ensure the central location of each 27.5 mm outer diameter latex cylinder 
within the associated 4-slit PVC pipe, each end of the latex cylinder was secured on to a Perspex 
cap which fitted inside the PVC pipe. Smaller PVC pipes with inner diameters of 22 mm were used 
inside the 4-slit pipes for some measurements. 

4.2 Measurement system and data analysis 

The sound source was a Bruel & Kjaer point source loudspeaker controlled by a Maximum-Length 
Sequence System Analyzer (MLSSA) system enabling determination of impulse responses. Mea- 
surements were taken of the insertion loss (IL) spectra for single cylinders and arrays of cylinders in 
an anechoic chamber. Figures [8|^ a) and (b) show an example measurement arrangement. Supports 
for the 2 m long cylinders were provided by holed wooden boards at the top and base of each array. 
The lattice constant (L) for the arrays of 4-slit cylinders and concentric cylinders was 80 mm. 

To maintain their shape and vertical orientation, the latex cylinders were slightly inflated above 
atmospheric pressure through a common pipe connecting to a small pump. During the array 
measurements, the receiver microphone was positioned 50 mm from the nearest face of the array 
but on the opposite side to the source. The loudspeaker was placed 1.5 m away from the array, 
such that the source-receiver axis was normal to the array orientation (see Figure El^a)). Both 
source and receiver were 1.2 m above the floor of the chamber which had sections removed to 
reduce unwanted reflections (see Figure EJ^b)). Insertion loss spectra were calculated by subtracting 
signals received without and with the cylinder array present but with the support structure in place 
on both occasions. 



5 Concluding remarks 

The acoustical properties of three forms of 2D resonators have been investigated: empty N-slit pipes, 
a concentric arrangement with rigid pipe inner and 4-slit pipe outer and a concentric arrangement 
of an elastic shell inner and 4-slit outer pipe. A theoretical formulation uses boundary conditions 
dependent on polar angle to represent arbitrary positioned slits. For simplicity, in the studied 
geometries the slits are positioned symmetrically. Jump boundary conditions imposed on the slit 
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interface are used to represent the solution inside the slits. The versatility of the proposed method 
can be applied to various type of 2D resonators with concentric multilayered/solid cylinder inner. 

It has been found theoretically and experimentally that increasing the number of slits in an 
empty pipe causes an increase in the frequency of the Helmhotz-type resonance. A low frequency 
approximation which models the slit cylinder by an equivalent fluid layer predicts that the frequency 
increase is proportional to the square root of the number of slits and this has been confirmed 
experimentally. The concentric arrangements result in resonances associated with both circular 
and annular cavities. With an inner elastic shell, an additional axisymmetric resonance of the shell 
is preserved but modified by the presence of the outer 4-slit pipe. Coupling between the components 
of the concentric arrangement results in shifts in the resonant frequencies corresponding to each 
element of the composite configuration. This is similar to the effect that is observed in mass- 
spring systems with multiple degrees of freedom. A low frequency approximation for the acoustical 
properties of the concentric arrangement with an elastic shell inner correctly predicts the observed 
frequency shift in the axisymmetric resonance. 

When used in periodic arrays the concentric arrangements with inner elastic shells and outer 
4-slit cylinders result in additional sound attenuation in the low-frequency range below the first 
Bragg band gap while still preserving the Bragg band gaps. This arrangement is more practical 
that the use of unprotected elastic shells so is potentially useful basis for a sonic crystal barrier 
design. 
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